Fermi liquid interactions and the superfluid density in d-wave superconductors 
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We construct a phenomenological superfluid Fermi liquid theory for a two-dimensional d-wave 
superconductor on a square lattice, and study the effect of quasiparticle interactions on the superfluid 
density. Using simple models for the dispersion and the Landau interaction function, we illustrate 
the deviation of these results from those for the isotropic superfluid. This allows us to reconcile 
the value and doping dependence of the superfluid density slope at low temperature obtained from 
penetration depth measurements, with photoemission data on nodal quasiparticles. 

PACS numbers: PACS numbers: 74.20.De,71.10.Ay,74.72.-h 



The high temperature superconductors appear to sup- 
port well-defined quasiparticle (QP) excitations at low 
temperatures (T rfni-Tc) as suggested by penetration 
depthtl, transportuEKa, and angle resolved photoemis- 
sion spectroscopy!] (ARPES) experiments. Low tem- 
perature superconducting (SC) state properties of the 
cuprates thus appear to be consistent with d-wave BCS 
theory with nodal QP excitations. However, the impor- 
tance of correlations, at low T is evident with under- 
doping: experiments^ show that the superfluid stiffness 
D s (T=0) ~x, and the QP weight; at (ir,0) diminishes on 
approaching the Mott insulatoiQ. In this paper we ad- 
dress the question of interaction corrections to the tem- 
perature dependence of D s (T) . 

The in-plane superfluid stiffness D S (T) = 
(c 2 d/4-7re 2 A 2 ), with d the mean interlayer spacing 
along the c-axis, can be directly obtained from mea- 
surements of the in-plane penetration depth, A(T\_ 
D S (T) is found to decrease linearly with temperatureu 
D S (T) = D s (0) - AT for T <g-TL, with a slope A which 
is nearly doping independenttru (or weakly decreasing 
but nonsingular) as x^O. 

Clearly the linear drop in D S (T) is due to thermally 
generated excitations which contribute to the normal 
fluid density. BCS theory with noninteracting QP exci- 
tations around the four d-wave nodes leads to the result 



D.{T)=DM-—^T. 

7T V 2 



(1) 



where vf is the Fermi velocity and vi is related to the 
slope of the SC gap via v 2 — (l/kF)dA(8)/<Mlg =n / i , at 
the nodal Fermi wavevector kp. Mesot et al.EEl obtained 
the nodal QP dispersion parameters vf and i>2 as a func- 
tion of doping from ARPES data on Bi2212, and com- 
pared the D s slope obtained from Eq.(^) with A measure- 
ments. They found that the slope estimated in this man- 
ner is too large by more than_a factor of two at optimal 
doping — the ARPES results^] of up =2.5 x 10 7 cm/sec 
and V2 = 1-25 x 10 6 cm/sec lead to an estimated slope 
dD s /dT = 0.77meV/K, while .the .slope obtained from 
penetration depth experiments^'Lali^ is approximately 
0.33meV/K. Furthermore, this discrepancy increases 



with underdoping since vi measured in ARPES decreases 
marginally leading to a slight increase in the estimated 
slope dD s /dT on underdoping, while the slope obtained 
from penetration depth experiments in Bi2212 decreases 
somewhat with underdopingo. This is in contrast to the 
rather striking agreement between estimates from ther- 
mal transport measurementaj and ARPESE3 for the ratio 
vf/v2~20 at optimal doping in Bi2212. 

Following Refs.li-3ii3, we attribute this discrepancy to 
residual QP interactions or Fermi liquid corrections. We 
use here a phenomenological superfluid Fermi liquid the- 
ory (SFLT) to explore the effects of lattice anisotropy 
on QP interactionsJn more detail than in earlier studies; 
(see, however ref£j). Some of the results obtained be- 
low were summarized without derivation in a conference 
reportO. 

We note that thermal phase fluctuation E are ignored 
here, since we have shown elsewhereO that a proper 
treatment of the long-range Coulomb interaction results 
in their contribution to D s (T) being subdominant to that 
of the nodal QPs. 

Superfluid Fermi liquid theory: Fermi liquid (FL) 
theory for a normal Fermi system is based on the exis- 
tence of well-defined (coherent) QP excitations which are 
adiabatic continuations of the single particle excitations 
of a free Fermi gas. While transport and ARPES ex- 
periments suggest that the normal state of optimal and 
underdoped high-T c SC's is not a FL, nevertheless, sharp 
QP peaks do appear all over the Fermi surface (FS) deep 
in the SC state (for T<^T C ). Naturally, one is then led 
to consider a description of the SC state and its low ly- 
ing QP excitations as an adiabatic continuation of a BCS 
state with Bogoliubov QP excitations. 

The approach advocated in refs.tSEj, and adopted be- 
low, assumes that such a SC state may be viewed as a cor- 
related FL in which a pairing interaction has been turned 
onEHl. In this case, one can use the superflui d Fe pai liq- 
uid theory (SFLT) developed many years agociEEa, and 
generalize it to the anisotropic case. 

For a normal Fermi system, the change in free energy 
due to a change in the QP momentum distribution Sn^ 



2 



takes the standard form 



Quasiparticle current renormalization: The QP 

energy £ k = ?k + Sk' /0**> k')<5"k', leads to the QP veloc- 



SF[Sn k ] = 2^ £££n k + - ^ /(k, k')<5n k <5n k , (2) ity v k = v° + £ k , V/(k, k')<5n k / . The total QP current 



k k,k' 

where £ k is the dispersion for the QP of momentum k in 
the absence of other QP's, /(k, k') is the Landau inter- 
action function, and Sn k = J2 a Sn k ^. We have ignored 
the spin-dependent part of /(k, k') in order to simplify 
the notation; the generalization with spin is straightfor- 
ward, but not relevant for the present discussion. Wc 
will refer to the QP's obtained by setting /(k, k') = 
in the above equation, as non-interacting QP's. The dis- 
persion for these QP's is £ k which does include the mass 
renormalization. 

We now use the above functional to calculate the su- 
perfluid stiffness at low temperatures, in two steps: (1) 
we calculate the diamagnetic response to a vector po- 
tential and (2) we calculate the renormalization of the 
current carried by the interacting QP's, relative to free 
QP's, and use this to compute the paramagnetic current 
correlator of the QP's. We next use the above quanti- 
ties as inputs to a Kubo formula in the QP basis, which 
allows us to determine the superfiuid stiffness D S (T). 

Diamagnetic term: Let n k be the unperturbed equi- 
librium QP distribution. In the presence of the vector 
potential n k — > "- k+eA / c leading to a shift of the momen- 



tum distribution <5n k 



diamagnetic termEl as the change 5F to order A 2 : 

+ i£/(k,k')A^vXv> k . (3) 



'k+eA/c 



We calculate the 



k,k' 



V^V^ denote derivatives with 



Here we set e = c = 1 

respect to k M and k^ respectively, where /j,,v = x,y and 
the sum over fi, v is implicit. The term linear in A van- 
ishes, since the integrand is odd in k, and we get 



SF = -A"A V 
2 



E'° 

k 



£k%,„n k +E^ k > k ') V M<V> k 



k,k' 



2 p 



(4) 



where K^v is the diamagnetic response. 

Given the jump discontinuity in the "normal" state, 
QP distribution at the FS, we use V M n k = -2u£ M <J(££), 
where the factor of 2 arises from summing over both 
spins. Using the definition v k = V£ k leads to 

k, v = 2 £ v^viMl)^ E /( k > ^'XAMimi ) 



k,k' 







(5) 



where K® v = 2 X]k u kV u k^(£k) 1S t ne diamagnetic term 



for non- interacting QP's. 



J is then ^k v k n k, which reduces to 

J = ]T v k 6n k - J2 /( k , k')£n k ,Vn k , 



(6) 



k,k' 



where we have used ^ k v k n k = in the first term, since 
the equilibrium QP population does not carry any cur- 
rent. In the second term, we have transferred the k- 
derivative from /(k, k') to n k , with n k w n k at this or- 
der. This relates the current carried by the interacting 
QP, to that carried by a non-interacting QP which only 
has a mass renormalization. 

To make further progress in the specific case of a d- 
wave SC, we note that the dominant excitations in the 
low temperature state are those near the gap nodes. We 
therefore restrict our attention to the renormalization of 
the current carried by the QP's at the 4 nodal points 
located at kf, with M = I.. A. Setting V p n k ~ 
— 2u kM 5(£ k ) as before, we find that the contribution to 
the current at the M-th node 



(7) 



l (M)Sn k (M) is the non-interacting 



J M (M) = J°(M) 

= J°(M)/3 

where J„(M) = v F 
QP current. In arriving at the above result, we have 
interchanged the k, k' labels in the second term, used the 
symmetry, /(k, k') = /(k',k), and there is no implicit 
sum over ^ in Eq. d?]). 

The superfiuid stiffness: From the Kubo formula, 
we find D>£ v = — A M „(q — > 0,iu) n = 0) where 

A M „(q,KJ„) = M (q, (-q, -«*>„)) is the current 

correlator and we take the transverse limit of q — > 0. 
In the QP basis, there are no excitations at T = 
and D^(T = 0) = K^ u . At low temperatures, there 
are nodal QP excitations and the current operator in 
A(q, iu; n ) has matrix elements between the ground state 
and these excited states. The current carried by the QP's 
is however renormalized by the factor j3 F which leads to 
A = /3?,A , with A being the correlator for the non- 
interacting QP's. The correlator A is easily evaluated 
within BCS theory using the dispersion £ k , and is linear 
in T at low temperature in a d-wave SC. Further, there 
are polarization effects by which the flowing QP's lead 
to an internal (fictitious) vector potential arising from 
the /(k, k'), in addition to the applied vector potentials. 
This effect is important close to T c when there are a large 
number of QP's, but it is unimportant at low tem pera. - 
ture when there are very few thermally excited QP'sPIr 3 !. 
The superfiuid stiffness in a d-wave SC at low T is thus 
given by 



D s {T)=a F K 



2\n2v F 



V-2 



T 



(8) 
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where K° — (1/d) TrK® in e?-dimensions, assuming cu- 
bic symmetry. We now proceed to discuss the FL correc- 
tions a F , (3 F in more detail. 

Isotropic limit: For an isotropic system v F and k F 
are independent of the location on the FS and to* = 
kp/vp is the effective mass. The Landau interaction 
/(k, k') = /(k ■ k') and depends only on the angle 
between the two momenta on the FS. Retaining only 
the single Landau parameter relevant for this discussion, 
/(k • k') = (dn / de)^ 1 Fi cos 8 , where cos6* = k • k' and 
(dn/de) — m* /it is the total "normal" state QP density 
of states for both spins. It is then easy to see that in 2D 



= 8^ — (I 

TO* 



■Ft/2) 



(9) 



where the 2D electron density n — k F /2ir. From Eq. (||), 
we thus find a F = (1 + F%/2). (For the special case of 
a Galilean-invariant system, using the Landau relation 
(l+Fi/2) = m*/m in 2D, we find K^ v = S^(n/m)). It is 
also easy to find that the renormalization of the current in 
the isotropic case is given by J p = J° (1 + i*\/2), and the 
current correlator is then A = 0* F A with (3 F = (l+Fi/2). 
These results for a F , [3 F ar&iri agreement with the earlier 
work of Larkin and MigdaEil and LeggettEa. 

We now discuss the shortcomings of isottopic SFLT as 
applied to the high T c SC's following Ref.E-3. Low tem- 
perature penetration depth experimentsQ suggest that 
D s (x,T = 0) ~ x. At the same time, ARPES experi- 
ments, asjstell as theoretical studies of SC in doped Mott 
insulatorsEj, suggest that m* does not diverge on un- 
derdoping. Within the isotropic SFLT framework, these 
two together imply (l + i<i/2)^x which in turn means 
the slope of D s (x,T) is -proportional to (1 + Fi/2) 2 ~x 2 . 
This scaling of the slopec3, however, is in strong disagree- 
ment with penetration depth measurements. Following 
the suggestions that this problem may be resolved by 
including anisotropy of the Landau interaction function 
over the FS, we next try to understand FL corrections in 
the anisotropic case. 

Anisotropic case: In order to set up a phenomeno- 
logical SFLT on a 2D square lattice, we first rewrite 
all our functions in terms of an angle variable 6 which 
sweeps over the large hole-barrel FS centered around 
(7r, 7r). Then, the Fermi momentum k F = k F (9), the 
Fermi velocity v F = v F (9) and the Landau interaction 
function /(k, k') = f(9,9'). We expand these in an or- 
thogonal basis 



v FX (0)=Y, V x )cos tt 2i+1 W 

oo 

v F A0) = J2 v r )s ' m tt 2i + 1 W 

£=0 

oo 

k F (9)=k FO +Y / ki e) co S (M9), 



(10) 

(11) 
(12) 



where we have used the symmetries of the square lat- 
tice to restrict the form of the expansion, and also 



used the vector (scalar) character of the v F (k F ). We 
may also generally expand the interaction, f(9,9') = 
Y,e,m F e,™e M e lme ' . We restrict the form of f{9,9') us- 
ing the following symmetries: (i) f(9,9') = f(9',9), (ii) 
f(0,6') - f{-9,~9') and (iii) /(tt/2 - 0,tt/2 - 9') = 
f(9,9'). While (i) is generally valid, (ii) and (iii) are 
valid for a square lattice. This finally leads to 



f(9, 6')=Y^ Ft,m [cos(£9 + me') + cos( 



m9)] (13) 



e>m 



where £,m : —oo — » oo with (£ + m) — 4jj and p — 
0, ±1, ±2, . . .. We have set £ > m to avoid overcounting. 
We note that: (a) the interaction function depends on 9 
and 9' separately in general and not only on (9 — 9') as 
in the isotropic case and (b) there are many more Lan- 
dau parameters on the lattice, labeled by two integers 
{£ , to). As we shall see, this considerably complicates our 
problem since many Landau parameters may contribute 
to a given response function, which prevents their unique 
determinationEH. This is unlike the isotropic case (say in 
He 3 ) where usually a single Landau parameter renormal- 
izes a particular correlation function. 

We now write the results for a F , f3 F in these new co- 
ordinates. The diamagnetic term is given by 



K<nx 2 



2tt 



d0 k F {9) 



(0) 



(14) 



. 2 *d9d9' k F {9) k F (9') 
o (2tt) 4 \v f (9)\ \v f (9')\ 



and the current renormalization for node-M is 
MM) [ 2 * d9> k F (9') v FX {0') m ^ n ,, 

JEm) + Jo W? \ v f W)\ v fx Wm ) [ ] 

where 9 M is the angular position of node-M. We can 
express this in a more compact form by defining (O)e = 
J 27I d9k F (9)0(9)/(2TT\v F (9)\). This yields 



a. 



I + ((v FX (9)v FX (9')f(9,9>)))ee> (lg) 



F = 1 



7r(v 2 )g 

(v F A0')f{0 M ,0'))e> 
™ FX (9 M ) 



(17) 



For f(9,9') = (7t/to*)Fi cos(9 - 9') and k F: v F indepen- 
dent of 9, we easily recover the isotropic limit. 

Simple models for the dispersion and /(k, k'): We 
now consider special cases of the general result which 
serve to illustrate the deviation from the isotropic limit. 

Case I: Consider an isotropic dispersion, with v F and 
k F independent of 9, but retain all allowed Landau pa- 
rameters on the lattice. In this case, with to* = k F /v F , 
we find 



a F = 1 + — (Fi,i + -Fx _!) 

7T 
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0.1 0.1 

Doping (x) Doping (x) 

FIG. 1: Doping dependence of the SFLT renormalization (f3 2 F ) 
of the slope of D S (T) for a model with anisotropic QP dis- 
persion and a single Landau parameter chosen such that (a) 
a F (x) = 1.5x and (b) a F (x) = 2.5a; (see Case II in text for 
details). In the isotropic limit, (3 F (x) = a 2 F (x), but there is 
marked deviation from this in the anisotropic case — most 
strikingly (3 F (x) t^D as x —* 0, as in the experiments. For 
this simple model and choice of dispersion, a larger renor- 
malization of Ds(0) (smaller a(x) as in panel (a)) appears to 
correlate with a weaker doping dependence of (3 F (x). 



\p<0 p>0 

+ £(-l)^4 P -i,i+£(-l) p F_ Mp+1 ) (18) 

p>0 p<0 / 

Thus, many Landau parameters contribute to the renor- 
malization in this anisotropic case unlike in the isotropic 
limit. Furthermore, different Landau parameters con- 
tribute to a F and (3 F . It is then easily possible that 



a F ^ [3 F and they could then also behave very differently 
with doping if several Landau parameters are nonzero. 

Case II: We next consider the case where we keep a sin- 
gle Landau parameter F\-\ ^ 0, and set all other Fg^ n — 
0. We however retain the full anisotropy of the disper- 
sion, as measured in ARPES. We take the_tight-binding 
fit to the (normal state) ARPES dispersioro, and numer- 
ically compute the above integrals to determine a F ,/3 F . 
In order to study the doping dependence of a F , (3 F , we 
assume a doping dependence Fi t -i(x) = B + Cx, such, 
that a F (x) ~ x in agreement with the Uemura plotfl, 
with reasonable values a F (x — 0.2) m 0.3 — 0.5. This 
fixes B, C and we use this to determine the doping depen- 
dence of f3 F (x) . The result of this calculation is plotted 
in Fig. |l|(a,b) where we see a marked deviation from the 
isotropic result ((3 F = a? ) — in the anisotropic case, f3 F 
is nonsingular as x — > 0, in qualitative agreement with 
penetration depth results. 

Conclusions: We have used a phenomenological SFLT 
for a d-wave SC to determine the renormalization of 
D S (T = 0) and dD s /dT due to FL factors. Within sim- 
ple models for the dispersion and the Landau interac- 
tion function, we find that anisotropy can cause strong 
deviations from the isotropic result. This allows us to 
understand the discrepancy between penetration depth 
and photoemission experiments for the temperature and 
doping dependence of the superfluid density in terms of 
SFLT corrections. While we discussed the case of a d- 
wave order parameter as appropriate for the high T c SC's, 
our results are easily generalized to any unconventional 
SC with point nodes and well-defined QP's. 
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